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We investigate the effects of quenched disorder on two chain Hubbard models at half-filling 
by using bosonization and renormalization group methods. It is found that the sufficiently strong 
forward scattering due to impurities and the random gauge field, which is generated by impurity 
backward scattering, destroy the charge gaps as well as the spin gaps. Random backward scattering 
due to impurities then drives the resulting massless phase to the Anderson localization phase. For 
intermediate strength of random forward scattering, however, the spin gaps still survive, and only 
one of the charge gaps is collapsed. In this parameter region, one of the charge degrees of freedom 
is in the Anderson localized state, while the other one is still in the massive state. 
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I. INTRODUCTION 

Coupled chain systems with ladder structure have at- 
tracted much current interest in the study of strongly 
correlated electron systems. It has been proposed that 
these systems QQirid show a superconducting state in the 
metallic phase.tl lI Such ladder systems have been reak 
ized experimentally in Srn-iCun+i02n and VO2P2O7.QEI 
More recently it has been reported that hole-doped 
Srn_iCun+i02n systcms actually undergo the transi- 
tion to the superconductivity under pressurcQ As mi- 
croscopic models for such systems, two coupled Hubbard 
chains cuad, Lwii, coupled t-J chains have been extensively 
studied.En3'Lrt3 The ground state of these models at half- 
filling is the Mott insulator with spin gap. When holes 
are doped into these systems, a metallic state with spin 
gap is realized. The resulting massless charge mode be- 
longs to a class of the Tomonaga-Luttinger (TL) liquid 
(sometimes referred to as the Luther-Emery class), and 
has the enhanced fluctuation toward superconductivity 
This superconducting state may be related with that 
found in Srn-iCun-i-i02n- 

Such ladder systems are quasi-one-dimensional, so that 
they may be quite sensitive to random potentials. It is 
thus important to investigate the effects of quenched dis- 
order on two coupled Hubbard chains. The case away 
from half-filling was studied by Orignac and Giamarchi 
beforeJlj However the case at half-filling has not been 
considered so far. In this case, we can systematically 
study how the Mott transition competes with the An- 
derson localization when the disorder is introduced to 
the system. This issue was previously addressed by the 
present authors, but for one-dimensional (ID) interacting 
electron systems.llZl(see also refs.18 and 19) In contrast 
to such single chain models, the Mott insulating state for 
two coupled Hubbard chains at half-filling has the spin 
gap as well as the charge gap. Since the randomness af- 
fects not only the charge gap but also the spin gap, it 



is expected that the presence of the spin gap may bring 
about novel properties, which have not been observed in 
the single chain model, for the competition between the 
Mott transition and the Anderson localization 

In this paper, we systematically study the effects of 
disorder on the two coupled Hubbard chains at half- 
filling. By exploiting bosonization and renormalization 
group methods, we discuss how the introduction of dis- 
order drives the system to the Anderson localized state. 
In particular, we focus on the role played by the spin gap 
for the competition between the Mott insulator and the 
Anderson localization. 

We deal with two types of quenched disorder: random 
impurity potentials within each band, and random hop- 
ping between two bands. We first discuss the effects of 
the former type of disorder. As was pointed out in ref. 17, 
the backward scattering due to impurities generates ran- 
dom gauge field coupled with electron currents. We find 
that the resulting random gauge fields play an essential 
role for destroying a charge gap in the coupled-chain sys- 
tems. At half-filling, sufficiently strong random forward 
scattering due to impurities and random gauge fields de- 
stroy charge gaps as well as spin gaps, and the resulting 
massless charge modes are localized by impurity back- 
ward scattering. However, for intermediate strength of 
random forward scattering, only one of the charge gaps 
is collapsed, and other modes are still gapful. In this 
case, we find a parameter region where one charge de- 
grees of freedom is in the Anderson localized state, while 
the other one is still in the massive state. 

In the case of random hopping between two bands, 
the analysis is more difficult, because two charge modes 
are mixed by random hopping. However introducing a 
unitary transformation, we can avoid this difficulty and 
find the same low-energy fixed point as the case without 
random interband hopping. 

The organization of the paper is as follows: In Sec. II, 
we introduce the model for the coupled Hubbard chains. 



1 



and apply bosonization methods to obtain low-energy ef- 
fective theory. For later convenience, we also summarize 
the results known for a clean system without random- 
ness. In Sec. Ill, we derive the scaling equations for a 
disordered system at half-filling. Sec. IV and V are de- 
voted to the discussions about the fixed point properties 
in the half-filling case. In particular, we discuss in detail 
how the competition between the Mott transition and the 
Anderson transition occurs when the random potentials 
are introduced. In Sec. VI, we further take into account 
the effects of random hopping between two bands. Con- 
clusion is given in Sec. VII. 



II. MODEL AND BOSONIZATION METHOD: 
SUMMARY OF THE RESULTS FOR CLEAN 
SYSTEMS 

We first summarize the results for a clean system with- 
out random potentials, and then investigate the effects of 
disorder on the resulting fixed points. Our model Hamil- 
tonian in the clean case is given by. 



H = -t'^ 4sqC*+isq + h.c. + U 



i.OL 



(1) 



where Cisa{c[sa) annihilation (creation) operator for 
electrons with spin s =11 and chain-index a = 1,2 at a 
site i, and riisa = cj^a'^isa- The last term is the hop- 
ping term between two chains. The kinetic energy part 
is diagonalized in terms of a new basis. 



Csljk) - Cs2{k) 
Csi{k) + Cs2{k) 

V2 



(2) 
(3) 



The dispersion relations in this basis are given by eik — 
—2tcosk + t± and £2fe = — 2tcosfc — t±, resulting in 
the two decoupled bands. For t± > 2t, there exists 
a band gap between these two bands. In the follow- 
ing we restrict our arguments to the interesting case of 
t± < 2t, where the chemical potential may cross the two 
bands both at and near half-filling. We apply abelian 
bosonizati| Oa rn cthods to this model to obtain the scaling 
equations.E30 In the presence of spin degrees of free- 
dom, non-abclian bosonization is genepaUj: more conve- 
nient for preserving SU(2) symmetry.r^lri However, in 
our model the backward scattering due to impurities is 
expressed in terms of the operators which are not in- 
cluded in the operator content of SU(2) Wess-Zumino- 
Witten model. Therefore we exploit conventional abelian 
bosonization methods by carefully dealing with SU(2) 
symmetry of spin degrees of freedom. Passing to the 
continuum limit, we linearize the dispersions around the 



Fermi momenta fcpi, 2 = | cos ^{{ztt± — n)/2t)\ {fi: chemi- 
cal potential), and divide the electron operators into left- 
and right-going parts, ipsi ^ exp{iy/4Tr(j)siL{x) + ikFix) + 
exp(— i\/47r0si/j(a;) — i/cpia^), etc. Defining the new boson 
phase fields. 
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(4) 
(5) 
(6) 
(7) 

with (j)si{2){x) = (t)si(2)L{x) + (i)si(2)R{x) {s =T,i), and 
their canonical conjugate momentum fields, H+p{x) — 
dxO+p{x) = dxi4>+pL - 4>+pr), etc., we write down the 
low-energy effective Hamiltonian, 



H — Hq + Hi + H2 , 



E 

a— ±,i/— p,cr 



rf2:[-^^((9a:0a^(a;))^ 



"-^KAx)]., 



dx 



a 



cos(\/87r0ip -I- Six) 



Ho = 



r dx 

+U / — cos(v87r02p + (^22;), 
J a 

— cos ■\/47r6'_p cos V^cb^a 
a 

+U2 / — cos Vivr^-n cos \/47r0+cr 
J a 

— cos \/47r0_p cos(V47r0+^ 
a 



(8) 



(9) 



(10) 



+ t/3 
+ C/4 

+ U7 



6x) 



dx 



cos V 4:n6-cr cos v incjj+a 



dx 
a 
dx 



/dx 
— cos Vivr^-cr cos(V47r(/)+p -I- Sx) 
a 

cos(\/47r0+p + Sx) cos \/47r0_cr 
cos \/47r0-cr cos V47r0+CT 

/dx 
— cos(\/47r0+p + Sx) cos V47r(/)+o- 
a 

/dx 
— COsV^9-aCOsV^0^p, (11) 
a 

where Si^2 — 4fcFi_2 ~ 2tt, and S — 2(fcFi + fcF2) — 2tt. 
Hi and H2 are, respectively, the umklapp scattering 
terms within each band and the other scattering terms of 
electron-electron interaction. For on-site Coulomb inter- 
action, initially Ui — U for i — 1 ^ 1 and Us — Ug — 0. 
We have dropped the terms which include an oscillat- 
ing factor exp(±2i(A:Fi — kF2)x), because these terms are 
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irrelevant in the long wave-length limit. We have also 
omitted the forward scattering terms of electron-electron 
interaction, which just renormalize the TL parameters. 

We first consider the case of half-filling 5 = 0. Then 
-ffi-term is irrelevant and can be discarded after the 
renormalization procedure because of the oscillating fac- 
tors. In general, the above Hamiltonian produces the 
mass gaps both for the (±)-charge sectors and (±)-spin 
sectors. In particular, it is to be noted that there are 
two possibilities for the gap formation of the (— )-spin 
sector; the spin gap is generated owing to locking of the 
or field. In this connection, we wish to mention 
that there is a special case for which some peculiar be- 
havior in the mass generation is observed, i.e. the case 
of Ui = U2 = Us = Ui = U5 = Ue = Ut = Us = Ug 
and -ftT-o- — Iri-Sfhere the model has special symmetry 
enhancement. Q'LJ This situation realizes in the case with 
a long-range electron-electron interaction between par- 
allel spins, H' = V J2i s='\i''^isni+is, where the initial 
values of Us and Ug are non-zero. To see the special 
symmetry, we rewrite the interaction term H2 as, 

H2 = Ui / — [(cos\/47r0+p -f cos V47r(?!)+cr 



cos 



cos 



cos 



cos 



'47r0_p) (cos V 47r0_cr + cos v ind-cr) 
^AttO^p cos V47r(f>+a- 
^47r0_p cos V47r(f>+p 

^<j)+„cosViTT(f>+p]. (12) 

Then the Hamiltonian is invariant under the dual trans- 
formation O-a- This is nothing but the Kramers- 
Wannier symmetry of the Ising model. Thus only the 
half part of the boson degrees of freedom in the (— )-spin 
sector is massive, and the other half-part is decoupled 
and forms a massless Pfpde of the Ising class with the 
central charge c = 1/2.QI1j The microscopic model given 
by eq.(|l|) does not have this special symmetry, so that 
we can proceed with the following analysis by assuming 
that all the low-energy excitations are gapful. 

We investigate the effects of disorder on these fixed 
points in the following sections. We note that the effects 
of quenched disorder on ID interactipg plpctron systems 
have been extensively studied before.cj Lj Our approach 
is a generalization of these previous studies. We intro- 
duce weak disorder potentials into the model (H), 



Hd;s = Hf+ Hb, (13) 



dx rj2ix)[^pl2L{x)^s2Lix) + tpl2R{x)i^s2R{x)] 



J dxY^ 773(0;) [f 



{i^llLi^)^s2Lix) 



5=n 

-'^l2Bi^)'^slBix)) + h.c], 



(14) 



Hb= fdxY, Mx)e^'''^'i^lLi^)AiR{x) + h.c] 

+ fdxJ2 Mx)e^'''''42Li^)^s2R{x) + h.c] 

s=n 



+i^l2Lix)^slRix)) + h.c] 



(15) 



where Hf is the forward scattering part, and iJf, is the 
backward scattering part. rii{x) and Ciix) (i = 1,2) 
are the random potential fields within each band for 
forward and backward scatterings, respectively. The 
?73 (x) and ^3 (x) terms represent the random hopping be- 
tween two splitted bands. We assume that these random 
fields obey the Gaussian distribution law, {r]i{x)r]j{x')) = 
Dn^SijS{x - x'), {^i{x)^*{x')) = D^.6ijS{x - x'), and 
D.^(^^i = -D,,(^)2- Wc first consider the case of Djjs = 
= in Sec. Ill, IV, and V, and take into account 
these random hopping terms in Sec. VI. 



III. SCALING EQUATIONS FOR DISORDERED 
SYSTEMS IN THE CASE WITHOUT RANDOM 
INTERBAND HOPPING 



In the absence of the 773 and ^3 terms, we bosonize the 
random potential terms (|l^ and (15), and average over 
the random fields using the replica trick as in ref. 27. 
Then the random potential terms in the action are given 
by, 



Sdi. 



j '^''j '^^ j dr'J2d.<l>\p{x,r)d^<l>ip{x,T')^^^ j '^'^ j '^^ j E ^-'^-p(^' 

j ^''j j dT'Y^^dr^+p{x,T)dr^\p{xy) 
~~~~ J y^'^ J dT''Y^-^dr(l)''_p{x,T)dr(l)Lp{x,T) 

y" da; y" dr y" dr' COS y/^ {x, t) +(j)'_„{x, t)) COS y/^{(t)\^{x,T) + (t)^_^{x,T)) 
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X COS y/^{^\p{x, t) + (jjtpix, t) - (jj^pix, r') - (t)^_p{x, t')) 
X COS ./^{(f>\p{x, t) - r) - t') + (fyLpix, t')), 



(16) 



where i, j are replica indices, and initially -D,,± = {D^i + 
D^2)/2. The Da± terms are generated in the. process of 
renormalization, though initially Da± = O.EZi The Da± 
term is equivalent to the interaction with random gauge 
fields, 



Hrg± = J dxA±{x)da:9±p, 



(17) 



where the random gauge fields A±{x) obey the Gaus- 
sian distribution law, {Aa{x)Af3{x')) = DAa5ap5{x — x') 
{a,P — ±). The random fields ri±{x) and A-d^^^cau be 
incorporated into the shift of the phase fields. tZiO Defin- 
ing the new fields, 



<p±p{x) = (j)±p{x) + r)±(x), 
V±{x) = ^^J- I dx'-q±{x'), 

V±p y TT J 
0±p{x) EE e±p{x)+A±{x), 

— # 

K±pV±p V TT 



A±ix) 



dx'A±{x'), 



(18) 
(19) 
(20) 
(21) 

where {ria{x)r]f3{x')) = Dna5ai35{x — x') {a, (3 — ±), we 
can absorb the rj± [x) and A± {x) terms into the Gaussian 
part of the Hamiltonian. Thus these terms do not affect 
the massless irLodes apart from a minor change in correla- 
tion functions.EIIEZI However, as we will see momentarily, 
these terms are important for suppressing the charge and 
spin gaps in the massive sectors. 

We now derive the scaling equations for the above ef- 
fective Hamiltonian. Since we have so many coupling 
constants, it is pretty difficult to analyze all scaling equa- 
tions even at one-loop level. We here ignore the renor- 
malizatiop_of velocities for simplicity. In some previous 
studiesElta, it was shown that even with this approxima- 
tion, the mass generation can be described rather well 
qualitatively for clean systems. This is partly because 
in the clean case the mass generation is mainly deter- 
mined by the scaling dimensions of the relevant operators 
which are not affected by the renormalization of veloci- 
ties in the lowest order. Thus this approximation gives 
a qualitatively correct answer for the mass generation 
in the clean case. We expect that this may be also the 
case so far as we are concerned with weak disorder. Us- 
ing the operator product expansion of the U(l) Gaussian 
model, and taking the replica limit n — > 0, we have the 
scaling equations for dimensionless couplings defined by 
Ui = Ui/v+p, D,j,A± = A,,A±/w±p, and D^i = D^i/v 



dUi 
~dl 



dU2 
dl 



dU^ 
dl 



^ = (2- 



dUi 
~dl 

dUs 
dl 



dUe 
dl 



dUj 

~dr 
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dl 
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dl 



dK. 



+p 



dl 
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"df 
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- K+„)Ui - 



K+p)U5 
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(2 - - K+p)Ue - 

2D^+KlpUe 
7r2 

(2 - - K+,)U7 ■ 
-D^i-D^2, 
(2 - - K+p)Us - 
2D,,+KlpUs 
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1 



K^p 



2DA-U9 



nKlp 



16 



(^51+^52), 



U2Ug 


U5US 
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2 




U3U9 
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U1U3 
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U1U2 
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2 
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U2U4 


U3U5 
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(22) 



(23) 



+p' 



16 



-iD^i + D^2), 



(24) 
(25) 

(26) 

(27) 
(28) 

(29) 

(30) 

(31) 
(32) 



4 



{Oil + ^52), 



dl 



16 



(33) 



dK-„ 7r,,-2 , , fV2l ^K-<y,fT2 



16 



dl), 



(3 



2 2 2 2 

UrD^„ 1 = 1,2. 



In the absence of disorder, _D^± = Da± = -D^i = 0, the 
excitation gaps open above the ground state energy in 
the (±)-charge and (±)-spin sectors. We have two pos- 
sibihties in the massive phase of the (— )-spin sector as 
mentioned in Sec. II. If Ui, Uq, and f/y scale to strong- 
coupling regime, and the Ui, C/5, and Ug terms are irrel- 
evant, the phase field (fi-a is locked, and the correlations 
for the operator exp(za0_o.) shows exponential decay. We 
refer to this case as case A. To the contrary, if the f/4, U5, 
and Ug terms are relevant, and the Ui, C/g, and f/y terms 
are irrelevant, is locked. We refer to this case as case 
B. In the subsequent sections, we discuss the effects of 
disorder on these two kinds of the fixed points. 



IV. DISORDERED FIXED POINT: CASE A 

Here we consider the effects of quenched disorder on 
the fixed point of case A, where the phase field c/)^^ is 
locked, and exp(ia9-a) is the disorder parameter. In this 
case, the C/i, U2, U3, U&, U7, and Us terms are relevant, 
and the U4, U5, Ug terms are irrelevant in the clean sys- 
tem. The sufficiently strong forward scattering due to 
impurities 1?,,+ may suppress the mass-generating inter- 
actions, C/3, Uq, and Us, as seen from Eqs.(p^, (|2^), and 
(p9|). It turns out from Eqs.(|2^) and ( p3| ) that this is not 
the case for the Ui and U2 terms which bear the charge 
gap in the (— )-sector, because they are not coupled with 
the random forward scattering. We note, however, that 
these terms should be suppressed in the presence of suffi- 
ciently strong random gauge fields A_ generated by ran- 
dom backward scattering due to impurities. Therefore in 
order to close the charge gap in the (— )-sector, one needs 
the random backward scattering due to impurities. This 
point makes a clear contrast to the case of a single chain 
model where the Mott-Hubbard gap is collapsed ouly by 
the random forward scattering due to impurities. tZI 

Note that the Ur term which generates the spin gap in 
(j)±a fields is not affected by the random forward scatter- 



ing, because the random forward scattering couples only 
with the charge degrees of freedom. On the other hand, 
it is seen from Eq. that the backward scattering due 
to impurities D^i may drive C/7 to a large negative value. 
However, if the initial values of D^i are small enough com- 
pared to that of Ur, the spin gaps of 0±o- fields generated 
by the Uj term may survive. Based on this observation, 
we start with the limiting case f/y ^ -foo for simplic- 
ity. Then after the spin gap formation, we can eliminate 
the spin degrees of freedom from the scaling equations. 
Suppressing irrelevant terms, we have the reduced scaling 



(34) 
(35) 
(36) 



(37) ^ ^ 



equations 
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U3U6 
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--K, 
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2DA-U2 
UiUe U5U9 



(38) 
(39) 
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2b^+KlpU3 



-± = {2-K+p)U^- 



dUe 
dl 

dU 
~dl 

dK+, 
dl 



dK_ 
"df 

db. 



{2~K+p)Us 



U1U3 
2 

2br, 



2Dr,+KlpU 



-KlpUs 



2 

ttK^ 



qu 

+p 



+ u' + u' 



16 



ib 



D 



«2J 



u' 



16 



dl 



= (3- 



i = 1,2. 



(40) 

(41) 
(42) 

(43) 
(44) 
(45) 



Since initially I?r)± > Da± = 0, and they obey the same 
scaling equations (35) and (^), flows into strong- 
coupling regime faster than Da±- Thus there exists a 
parameter region where the C/3, Uq, and Us terms are 
suppressed for sufficiently large , though the Ui and 
U2 terms are still relevant. Then the charge gap in (-f )- 
charge mode is collapsed. Since the 9-p field is locked 
due to the Ui and U2 terms, K^p scales to a large value, 
and the b^i terms are suppressed, as seen from Eq.(^5l). 
However, these terms are proved to be relevant, after one 
integrates out the (p-p ifipld, following the argument by 
Orignac and GiamarchiE3. as in the case away from half- 
flling£3 Then the (-l-)-charge mode is in the Anderson 
localized state, while the (— )-charge mode is still in the 
Mott insulating state. 

As the value of bA- becomes larger, the Ui, U2 terms 
are suppressed, and all excitation gaps are collapsed. 
Then, the backward scattering due to impurities, 
is relevant, and consequently the Anderson localization 
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occurs. Ur flows to a large negative value because of the 
D^i term, as seen from Eg. (Pq) . Thus the spin degrees 
of freedom are frozen. The fixed point is identifled with 
the pinned CDW state which is sijiailar to that found in 
the disordered single chain model J23 

In order to confirm the above arguments, we solve the 
scaling equations (^2|)-(^7|) numerically for several initial 
values of coupling constants. The numerical results for 
the flow of coupling constants are shown in FIG.l. In 
order to realize the massive spin state due to the lock- 
ing of in the clean system, we have chosen the ini- 
tial value of the TL liquid parameter in the (— )-spin 
sector as K-o-{0) = 0.86. The other TL parameters 
are set to be K±p{0) = 0.98 and K+„{0) = 1.0. For 
these initial values, the system definitely flows toward 
the low-energy fixed point of case A in the absence of 
disorder. In FIG. 1(a) the results in the case without 
quenched disorder are shown. Ui, U2, f/3, C/g, Ur, and 
Us flow into a strong-coupling regime, while C/4, C/5, and 
C/g scale to 0. As a consequence, the phase fields (/)+p, 
8^ p. (j)±cr are locked. This fixed point is stable against 
the disorder weaker than mass-generating interactions, 
because of the presence of excitation gaps in all modes. 
As the quenched disorder becomes stronger, some differ- 
ent fixed points manifest themselves. In FIG. 1(b), we 
show the results for the initial values I?,,±(0) — 0.9 and 
D^i{0) = 0.005. Although Ui, U2 and U7 still scale to 
a strong-coupling regime, U3, Uq, and Us scale to 0. As 
a result, the charge gap in the (-l-)-mode is collapsed. 
Since D^^i is relevant, the charge mode in the (-l-)-sector 
is localized, though the charge mode in the (— )-sector 
is still massive. Thus the Anderson localized state and 
the Mott insulating state coexist in this phase. Since 
C/7 —f —00, the spin degrees of freedom are frozen, and 
may be in the random spin singlet state. In FIG. 1(c), the 
results for Djj±{0) = 0.9 and -D{i(0) = 0.01 are shown. 
In this case, all mass-generating interactions except U-j in 
Eq.(pTl) flow to 0. The charge gaps are thus collapsed in 
this phase. Since D^i scales to a strong-coupling regime, 
the Anderson localized state realizes. The relevance of 
the t/7 term implies freezing of spin degrees of freedom 
in the (±)-sectors. 

In summary, we have found three kinds of fixed points 
in the presence of quenched disorder: (i) the clean fixed 
point where the disorder effect is entirely irrelevant, (ii) 
Anderson-insulator (I) where the (-l-)-charge mode is in 
the Anderson localized state, the (— )-charge mode is in 
the Mott insulating state, and the spin degrees of freedom 
are massive or frozen, (iii) Anderson-insulator (II) where 
all charge modes are in the Anderson localized state, and 
spin degrees of freedom are frozen, i.e. pinned CDW. We 
thus end up with the schematic phase diagram for charge 
degrees of freedom as shown in FIG. 2. 



V. DISORDERED FIXED POINT: CASE B 

Here we consider the effects of disorder on the fixed 
point of case B: d^a is locked, and exp(ja(/)_(j) is the 
disorder parameter field in the absence of quenched dis- 
order. In this case, the [/2, t/s, C/4, C/5, C/s, and C/g terms 
are relevant in the clean system. We see from Eq.(|3^) 
that K-a is decreased by the backward scattering due 
to impurities Z?^i. Thus it suppresses the growth of C/4, 
C/5, and C/g and prevents the spin gap formation due to 
locking of the 6-^ field. Moreover it follows from Eq.([2^) 
that if the D^i term is relevant, C/7 flows to a large neg- 
ative value, resulting in freezing of the spin degrees of 
freedom, and the suppression of the C/4 term as in Sec. 
IV. However for sufficiently small D^i, the C/4 term is rel- 
evant and the spin gap opens owing to locking of 0-„. 
After the generation of the spin gap due to locking of the 
6'_o- and 4>+cr fields, we can eliminate the spin degrees of 
freedom, and obtain the scaling equations for the charge 
degrees of freedom. 



dU_2 
dl 

dU3 
dl 



1 



)C/2 



U3US 



(46) 



2Da^U3 2D^+Kl„U3 



^ = (2-i.,,)C/5 



U3U9 



2Dr,+ KlpU5 



dUs 
dl 

dUg 
dl 



2D^+KlpUs 



(2 - K+p)Us - 2 

TT 



(2- 



1 ,~ U3U5 2DA-Ug 

1^9 



dK^ 



+p 



dl 



dK^r 



ttKIp 
16 



'-[ui + u^ + uil 

{D^i+D^2), 



dl 



;[U2+U! + U^]. 



'-{D^,+D^2), 



dDn 



dl 



16 

(3 - 2K+P - 2K.p)b^, 



(47) 

(48) 
(49) 
(50) 

(51) 

(52) 
(53) 



As pointed out in Sec. IV, Dr)± grows faster than Da±- 
Thus there may exist a parameter region where D,,+ sup- 
presses U3, C/5, and C/g, but the U2 and C/g terms are still 
relevant, bearing the mass gap in the (— )-charge mode. 
In this parameter region, only the charge gap in the (-1-)- 
charge sector is collapsed. The low-energy properties of 
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this sector is determined by the D^i-terms. After inte- 
grating out the massive (— )-charge modell3, we find that 
the scaling equation for the -D{i,2-terms are changed to, 



dD, 



dl 



(3 - 8K+p)D^, 



(54) 



Thus the impurity backward scattering is suppressed in 
comparison with single chain systems. For if+p > 3/8, 
the Z3^i_2-terms are irrelevant. However it does not di- 
rectly mean the occurrence of delocalization, since we 
should actually include the random hopping between two 
bands, D^3, D^^. In the case of D^a ^ and D^^ ^ 0, 
the CDW order couples directly to the D^^ term, and 
then the pinning of the CDW occurs for infinitesimally 
small Z?^3 as will be discussed in Sec. VI. Therefore the 
(-l-)-charge mode is in the Anderson localized state with 
the vanishing Drude weight and the non-zero charge sus- 
ceptibility. For a sufficiently large value of Da-, the 
charge gap in the (— )-mode may be also collapsed. All 
the charge modes are in the Anderson localized state. 

In order to specify the low-energy fixed points cor- 
rectly, we have solved the full scaling equations (|2^)-(|37|) 
numerically for some initial values of couplings. The nu- 
merical results are shown in FIG. 3. In FIG. 3(a), we 
display the results for clean systems. In order to re- 
alize the spin gapped state due to locking of 0-a, we 
choose the initial values of the parameters K-cr{0) = 1-1, 
K-^-cr{0) — 1.0, a.nd K±p{Q) — 0.98. For these parameters, 
U2, U3, C/4, f/5. Us, and f/g scale to a strong-coupling 
regime, and t/i, Uq, and scale to 0. Then the phase 
fields (p+p, 0-p, (f)+a, and 6-a- are locked. This fixed 
point with the spin and charge gaps is stable as far as 
disorder is much weaker than mass-generating interac- 
tions. In FIG.3(b), the results for i?^±(0) = 0.9, and 
D^i{0) — 0.007 are shown. This flow ends up at the 
same fixed point as that shown in FIG. 1(b), where the 
backward scattering due to impurities is relevant, and 
thus the Anderson localization occurs in the (-l-)-charge 
mode, though the (— )-charge mode is in the Mott insu- 
lating state. The spin degrees of freedom are also massive 
or frozen. Note that the freezing of the spin degrees of 
freedom in the (— )-sector is not due to locking of 0„o- 
but due to (/"-o-. Thus the fixed point properties of the 
(— )-spin sector are different from those for the clean case 
where locking of d-a occurs. In FIG. 3(c), we show the 
results for D,i±{0) — 0.9 and D^i{0) — 0.01. In this case, 
the flows of coupling constants are qualitatively the same 
as those shown in FIG. 1(c); i.e. all excitation gaps are 
collapsed, and the Anderson localized state with frozen 
spin degrees of freedom is realized. Thus for sufficiently 
strong impurity scattering, the disordered fixed point of 
case B flows into the same class as that of case A. In 
summary, we have found three kinds of fixed points as 
in case A. The schematic phase diagram for the charge 
degrees of freedom is also given by FIG. 2. 



VI. EFFECTS OF RANDOM HOPPING 
BETWEEN TWO BANDS 

We now take into account the effects of random hop- 
ping between two bands, 773 and ^3. First we consider 
only random forward scattering, and omit random back- 
ward scattering for a while. Going back to the ferniion 
representation, we write down the kinetic energy part and 
the impurity forward scattering part of the Hamiltonian, 

s=Ti 

+ V2{iP12l9x^s2L - 1pl2Rdxi^s2R)] 

dx r]2ix)[ijl2L{x)^,2Lix) + i^l2ji{x)i:s2Rix)] 

+ fdxY, m{xW^'-^-''-^>{i^l^{x)i;s2L{x) 

+'^l2Ri^)'^slR.{x)) + h.c.]. (55) 

In the case of half-filling, vi — V2 ^ v holds. Defining 
the spinor fields, 



'^sL{R) 



i'slL(R) 



(56) 



'0s2L(i?) 

we cast the Hamiltonian Eq. ( |55| ) into the following form. 



Hq = J dx Y vitj^l^dxil^sL - i^lRdx^JsR) 



s=n 



1 H T 



/ dxY^^\^[ 



V1+V2-, , V1-V2 z 

:^ 1 H :^ T 



where 1 is a 2 x 2 unit matrix, and r^, are the Pauli 
matrices. This Hamiltonian can be diagonalized by using 
the unitary transformations defined by, 



IpsL = UliPsL, 

Ul = expf / dx' 

, r]i{x') - mix') 



(58) 



rii{x') + 772(2;') . 



-f?73(a;')e*^'=^i-'=^^)^' — 
+77*(2;')e-*('="-'=''^)^'^ ) 



(59) 
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and Ur = Uj^. We thus end up with the diagonahzed 
form of Hq, 



VII. CONCLUSION 



(60) 



Therefore the forward scattering part of the random hop- 
ping can be incorporated into the free part of the Hamil- 
tonian. Under the unitary transformation (^9|), the for- 
ward scattering terms of electron-electron interaction are 
invariant. The backward and umklapp scattering terms 
of electron-electron interaction ( [Til ) which may generate 
the mass gaps are multiplied with exponentially decay- 
ing factors by this transformation. Thus for sufficiently 
strong impurity forward scattering both of the charge 
and spin gaps collapse as in the previous sections. 

Let us now take into account random backward scat- 
tering. The backward scattering term of random hopping 
is expressed as 



Hi3 - / da;[e3(a;)e^('=- 



OcDw{x) + h.c]. (61) 



The scaling equation for D^^ is given by 

1 K+, 1 



dl 2 



2K- 



2K^ 



(62) 



where D^^ = D(^^/v\p. It is seen that the D^^ term cou- 
ples with the J74 term. Then the scaling equation (|2~' 
altered to 



IS 



dUi 
"df 



1 



C/4 - U2U9 - U5US 



-D. 



(63) 

Therefore in the case that the mass gaps open due to 
locking of the 0+p, 0-p, and fields in the clean 
system (case 5), ■i-K+p/2-l/2K-p-K+a/2-l/2K-a- 
[/4 > holds, and the term is relevant. Then, the 
Anderson localization inevitably occurs in the (-l-)-charge 
sector provided that the impurity forward scattering sup- 
presses the charge gap as mentioned in Sec. V. This fixed 
point is the pinned CDW. If the random backward scat- 
tering D^i 2 drives C/7 to —00, the 0_cr field is locked. 
Then in the presence of the gap in the (— )-charge sector, 
integrating out the 6^?ie\d as in ref. 16, we find that 
the scaling equation (62) is changed to. 



^ = (3-2i^+,)i?,3. 



(64) 



Thus although the Z3^3-term is suppressed, it is still rel- 
evant for the repulsive interaction Kj^p < 1, and the An- 
derson localization occurs in the (-l-)-charge sector. 

If the spin gap opens owing to locking of the 0-^ field 
in the clean system (case A), the scaling equation for 
-Dj3 is also given by Eq.(^J). Then the random hopping 
between two bands enhances the Anderson localization 
obtained in Sec. IV. 



We have systematically studied the effects of quenched 
disorder on the spin-gapped state of two chain coupled 
Hubbard ladder model at half-filling. We have shown 
that a sufficiently strong forward scattering due to impu- 
rities and the random gauge fields generated by impurity 
backward scattering collapse all the charge gaps as well 
as the spin gaps, and consequently the Anderson local- 
ization takes place in all the charge sectors owing to the 
impurity backward scattering. We have also found the 
phase that one charge mode is in the Anderson localized 
state, while the other one is in the massive state. 

This work was partly supported by a Grant-in-Aid 
from the Ministry of Education, Science and Culture. 
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FIG 1. Plots of Ui and vs the scaling parameter 
I in the case A for different initial values of Djj± and 
Di^■■ (a) -D^±(0) = D^,iO) = 0. (b) -D^±(0) = 0.9, 
D^i{0) = 0.005. (c) A,±(0) = 0.9, D^i{0) = 0.01. The 
initial values of the TL liquid parameters and Ui are set 
as K±p{0) = 0.98, ^+^(0) = 1.0, K_„{0) = 0.86, and 
Ui{0) = 0.1 for i = 1 - 7, jJsiO) = UgiO) = 0. As the 
initial values of and D^i increases, the mass generat- 
ing terms scale to zero, and the Anderson localized state 
realizes. 
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FIG 2. Phase diagram in the Dri±-D^i plane, (i) Mott 
insulator (MI) in both the (+) and (— )-charge sectors, 
(ii) Anderson localization (AL) in the (-l-)-charge sector, 
and MI in the (— )-charge sector, (iii) AL in both the (+) 
and (— )-charge sectors. 
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FIG 3. Plots of Ui and Z?^ vs I in the case B. for the 
different initial values of I?^± and .0^^: (a) -D^±(0) = 
^^,(0) = 0. (b) ^^±(0) = 0.9, ^4,(0) = 0.007. (c) 
i)^±(0) = 0.9, i>^i{^) = 0.01. The initial values of 
the TL parameters and Ui are set as i^±p(0) = 0.98, 
i4:+<,(0) = 1.0, i^-<,(0) = 1.1, {/j(0) = 0.1 for i = 1 - 7, 
and i/8(0) = 1/9(0) = 0. 
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